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The  Effects  of  Spatial  Diversity  and 
Imperfect  Channel  Estimation  on 
Wideband  MC-DS-CDMA  and  MC-CDMA 

Andrew  S.  Ling,  Member,  IEEE,  and  Laurence  B.  Milstein,  Eellow,  IEEE 


Abstract — In  our  previous  work,  we  compared  the  theoretical 
hit  error  rates  of  multi-carrier  direct  sequence  code  division 
multiple  access  (MC-DS-CDMA)  and  multi-carrier  code  division 
multiple  access  (MC-CDMA).  To  ensure  a  fair  comparison,  we 
constrained  both  schemes  to  the  same  bandwidth,  information 
rate,  and  energy-per-bit,  and  these  constraints  resulted  in  a 
possible  trade-off  between  diversity  gain  and  channel  estimation 
errors  between  the  two  schemes.  While  only  a  single-input  single¬ 
output  (SISO)  system  was  considered  in  our  previous  work, 
in  this  paper,  we  extend  the  comparison  to  a  multiple-input 
multiple-output  (MIMO)  system  which  employs  Alamouti  space- 
time  block  coding  at  each  sub-carrier  frequency  to  achieve  spatial 
diversity.  We  consider  only  those  cases  where  MC-CDMA  has 
higher  frequency  diversity  than  MC-DS-CDMA.  Since  increases 
in  diversity  yield  diminishing  gains,  we  conclude  that  the  addition 
of  spatial  diversity  to  this  multi-carrier  comparison  benefits  MC- 
DS-CDMA  more  than  MC-CDMA.  To  determine  whether  these 
gains  for  MC-DS-CDMA  are  enough  to  offset  the  difference  in 
frequency  diversity  between  the  two  schemes,  we  derive  closed- 
form  expressions  for  the  bit  error  probabilities  of  both  schemes, 
and  we  compare  the  MIMO  results  against  those  of  the  SISO 
system  for  different  information  rates,  number  of  users,  and 
number  of  pilot  symbols  per  channel  estimate. 

Index  Terms — MC-DS-CDMA,  MC-CDMA,  Alamouti  space- 
time  block  coding,  frequency-selective  Rayleigh  fading,  Herml- 
tian  quadratic  form. 


1.  Introduction 

MULTI-CARRIER  signaling  has  been  an  active  area 
of  research  over  the  past  fifteen  years.  Two  multi¬ 
carrier  schemes  which  have  been  widely  studied  are  multi¬ 
carrier  direct  sequence  code  division  multiple  access  (MC- 
DS-CDMA)  [l]-[7]  and  multi-carrier  code  division  multiple 
access  (MC-CDMA)  [8]-[12].  When  viewed  in  the  frequency 
domain,  these  two  schemes  differ  in  the  widths  of  their  sub¬ 
bands:  MC-DS-CDMA  uses  direct  sequence  spreading  at  each 
sub-carrier,  while  each  sub-carrier  in  MC-CDMA  is  unspread. 
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As  a  result,  over  a  given  bandwidth,  MC-CDMA  employs 
a  larger  number  of  sub-carriers  than  MC-DS-CDMA.  On 
the  one  hand,  this  implies  that  MC-CDMA  potentially  has 
higher  frequency  diversity  than  MC-DS-CDMA  when  both 
schemes  transmit  at  the  same  information  rate;  on  the  other 
hand,  this  also  implies  that  the  energy  per  sub-carrier  is  lower 
in  MC-CDMA  than  in  MC-DS-CDMA  when  both  schemes 
use  the  same  energy-per-bit.  Since  each  MC-CDMA  sub¬ 
carrier  operates  at  a  lower  signal-to-noise  ratio  (SNR),  the 
receiver’s  estimate  of  the  channel  gain  at  each  sub-carrier 
frequency  in  MC-CDMA  is  more  prone  to  error.  Thus,  there 
exists  a  possible  trade-off  between  diversity  gain  and  channel 
estimation  errors  between  these  two  schemes.  There  are  many 
papers  which  examine  this  trade-off,  as  well  as  other  aspects, 
when  comparing  MC-DS-CDMA  and  MC-CDMA.  To  the  best 
of  our  knowledge,  however,  the  numerical  results  of  these 
papers  are  all  based  on  computer  simulations  [13]-[19]. 

In  [20],  we  specifically  considered  the  problem  of  whether 
a  multi-carrier  system  should  use  spread  or  non-spread  sub¬ 
carriers.  Our  approach  to  this  problem  involved  a  comparison 
of  the  theoretical  bit  error  rates  between  MC-DS-CDMA  and 
MC-CDMA  and  took  into  account  the  aforementioned  trade¬ 
off  between  diversity  gain  and  channel  estimation  errors.  We 
were  not  concerned  with  optimizing  performance  or  adhering 
to  conventional  implementations.'  Rather,  our  main  objective 
was  to  establish  a  fair  comparison.  Therefore,  we  required 
that  both  schemes  use  (1)  the  same  bandwidth,  information 
rate,  and  energy-per-bit;  (2)  the  same  combining  scheme  at 
the  receiver  (i.e.,  maximal  ratio  combining);  and  (3)  the  same 
pilot-based  scheme  to  estimate  the  channel  gains  at  each  sub¬ 
carrier  frequency.^  We  then  derived  closed-form  expressions 
for  the  bit  error  probabilities  of  both  schemes  for  two  separate 
channel  scenarios  (i.e.,  two  different  cases  for  the  coherence 
bandwidth  of  the  channel),  and  we  compared  these  error 
probabilities  to  determine  the  performance  trade-offs  that  may 
result  from  using  one  scheme  over  the  other. 

In  this  paper,  we  extend  our  previous  comparison  between 
MC-DS-CDMA  and  MC-CDMA  to  a  multiple-input  multiple- 

*Our  MC-DS-CDMA  system  was  based  on  the  one  proposed  in  [2]  and 
studied  further  in  [3]-[5].  Our  MC-CDMA  system  was  modeled  after  our  MC- 
DS-CDMA  system  but  still  kept  the  main  essence  of  MC-CDMA — namely, 
that  the  data  at  each  MC-CDMA  sub-carrier  is  unspread  and  that  each  sub¬ 
carrier  is  multiplied  by  a  different  chip  in  the  spreading  sequence. 

^For  both  multi-carrier  schemes,  we  assumed  that  the  receiver  has  no 
knowledge  of  the  coherence  bandwidth  of  the  channel  and  is  forced  to  estimate 
the  channel  gains  at  each  sub-carrier  frequency. 
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Fig.  1.  Transmitter  block  diagram. 


output  (MIMO)  system.  For  both  multi-carrier  schemes,  we 
use  Alamouti  space-time  block  coding  at  each  sub-carrier 
frequency  to  achieve  spatial  diversity  in  the  absence  of  channel 
state  information  (CSI)  at  the  transmitter  [21].  By  assuming 
(1)  the  coherence  bandwidth  of  the  channel  is  equal  to  the 
bandwidth  of  one  MC-DS-CDMA  sub-band,  and  (2)  the 
number  of  data  symbols  transmitted  in  parallel,  M,  is  such  that 
M  >  1,  we  consider  only  those  cases  where  MC-CDMA  has 
higher  frequency  diversity  than  MC-DS-CDMA,  because  these 
are  the  only  instances  in  which  there  is  a  performance  trade¬ 
off  between  the  two  schemes  [20].  Since  increases  in  diversity 
yield  diminishing  gains,  we  conclude  that  the  addition  of 
spatial  diversity  to  this  multi-carrier  comparison  benefits  MC- 
DS-CDMA  more  than  MC-CDMA.  To  determine  whether 
these  gains  for  MC-DS-CDMA  are  enough  to  offset  the 
difference  in  frequency  diversity  between  the  two  schemes,  we 
use  a  quadratic-form-based  technique  to  derive  closed-form 
expressions  for  the  bit  error  probabilities  of  both  schemes, 
and  we  compare  the  numerical  results  of  the  MIMO  system 
against  those  of  the  SISO  system  for  different  information 
rates,  number  of  users,  and  number  of  pilot  symbols  per 
channel  estimate. 

This  paper  is  organized  as  follows.  In  Section  II,  we  present 
the  system  model  assuming  a  general  multi-carrier  signaling 
scheme.  We  apply  MC-CDMA  and  MC-DS-CDMA  signaling 
to  this  system  model  in  Section  III,  and  we  analyze  their  error 
probability  performances  in  Section  IV.  The  numerical  results 
of  the  MIMO  system  are  then  compared  against  those  of  the 
SISO  system  in  Section  V,  and  Section  VI  concludes  the  paper. 

Notation:  Lowercase  (uppercase)  boldface  letters  denote 
column  vectors  (matrices);  I^r  and  Oat  are  the  N  x  N 
identity  and  zero  matrices,  respectively;  (•)^,  (•)*,  and  (•)^ 
denote  transpose,  complex  conjugate,  and  complex  conjugate 
transpose,  respectively;  0  denotes  the  Kronecker  product; 


Sn,m  =  5[n  —  m],  where  <)[•]  is  the  Kronecker  delta  function. 

II.  General  MIMO  Multi-Carrier  System  Model 
A.  Overview 

Consider  an  uplink  scenario  involving  K  asynchronous 
users.  Assume  two  antennas  at  each  transmitter  and  two  an¬ 
tennas  at  the  receiver.  The  available  bandwidth  is  the  same  for 
each  user  and  is  divided  into  S  bandlimited,  non-overlapping 
frequency  sub-bands  of  equal  width. ^  We  assume  each  sub¬ 
band  is  associated  with  a  sub-carrier  frequency  and  has  a 
bandwidth  of  W . 

In  the  time  domain,  the  transmissions  made  at  each  sub¬ 
carrier  are  structured  into  frames,  each  consisting  of  Q  time 
slots.  As  in  [20],  we  assume  a  slowly-varying  channel  which 
remains  constant  over  the  duration  of  a  frame,  and  we  desig¬ 
nate  the  first  Qp  time  slots  of  each  frame  as  the  estimation 
phase  and  the  remaining  Qp,  =  Q  —  Qp  time  slots  as  the  data 
phase.  During  each  time  slot  of  the  data  phase,  both  antennas 
simultaneously  transmit  M  distinct  binary  data  symbols  in 
parallel  across  the  S  sub-bands,  where  each  symbol  is  repeated 
on  R  sub-carriers.  Thus,  we  have  S  =  MR,  and  the  overall 
system  bandwidth  is  equal  to  SW  =  MRW.  Similarly,  during 
the  estimation  phase,  we  transmit  M  pilot  symbols  (-fl’s)  in 
parallel,  and  we  spread  the  energy  of  each  pilot  symbol  across 
R  sub-bands.  However,  pilot  symbols  are  transmitted  on  the 
first  antenna  only  during  the  odd  time  slots  and  on  the  second 
antenna  only  during  the  even  time  slots.  As  a  result,  only 
the  odd  (even)  time  slots  are  used  to  estimate  the  channel 
gains  associated  with  the  first  (second)  transmit  antenna,  and 
each  estimate  is  formed  from  Qpj2  pilot  symbols."^  To  keep 

^Spectral  efficiency  is  not  a  concern  in  this  comparison,  so  we  do  not  insist 
on  using  minimum  sub-carrier  spacing  and  having  the  sub-bands  overlap. 

^We  assume  Qp  is  even. 
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the  total  energy  transmitted  from  the  two  antennas  the  same 
over  all  of  the  time  slots,  we  also  use  a  different  transmit 
amplitude  during  the  estimation  phase.  Our  motivation  for 
imposing  this  energy  constraint,  as  well  as  for  using  the  on-off 
pilot  transmission  scheme,  is  to  maintain  the  same  energy-per- 
pilot  symbol  between  this  system  and  the  SISO  system  in  [20], 
By  doing  so,  channel  estimates  formed  from  the  same  number 
of  pilot  symbols  will  be  of  the  same  quality  in  both  the  SISO 
and  MIMO  systems. 

B.  Transmitter 

The  transmitter  block  diagram  is  shown  in  Fig.  1.  During 
the  estimation  phase,  the  on-off  transmission  of  pilot  symbols 
on  all  the  sub-carriers  is  modeled  using  the  switches  at  the  two 
antennas.  The  switch  for  the  first  (second)  antenna  is  closed 
during  the  odd  (even)  time  slots  and  open  during  the  even 
(odd)  times  slots.  The  details  of  the  multi-carrier  modulation 
blocks  will  be  given  in  Section  III  when  we  consider  MC- 
CDMA  and  MC-DS-CDMA. 

During  the  data  phase,  the  input  data  is  modeled  as  an  inde¬ 
pendent  random  binary  sequence  with  bit  rate  1/T.  The  serial- 
to-parallel  (S/P)  converter  produces  M  parallel  sequences,  and 
we  label  each  group  of  M  symbols  as  where 

p  is  the  frame  index,  q  is  the  time  slot  index,  and  k  is  the 
user  index.  A  rate  1/R  repetition  code  is  applied  to  each  of 
these  M  symbols,  and  the  resulting  MR  symbols  are  mapped 
to  the  MR  sub-carriers  in  such  a  way  that  maximizes  the 
frequency  separation  between  adjacent  repetitions  of  a  symbol. 
We  denote  the  frequency  to  which  the  i-th  repetition  of  the 
TO-th  symbol  is  mapped  by  After  symbol  mapping, 

Alamouti  space-time  block  coding  is  performed  at  each  sub¬ 
carrier  frequency,  and  the  outputs  of  the  encoders  are  fed 
through  the  two  banks  of  modulators.  Finally,  the  signals  are 
combined  and  transmitted  over  the  two  antennas. 

To  represent  the  transmitted  signals  at  each  frequency 
during  the  estimation  and  data  phases  under  a  common  mathe¬ 
matical  framework,  we  define  b^p}q^u,m  as  the  fc-th  user’s  input 
to  the  modulators  at  frequencies  {fm,i}i=i,...,R  of  the  rx-th 
transmit  antenna  during  the  g-th  time  slot  of  the  p-th  frame. 
As  shown  in  Table  I,  for  q  =  1, . . .  ,Qp,  is  equal  to 

either  -fl  or  0,  while  for  q  =  Qp  +  1,  ■ .  ■  ,Q,  is  equal 

to  the  output  of  the  Alamouti  encoder.  The  transmitted  signal 
for  the  /c-th  user  is  given  by  the  2x1  vector 

M  R 

=  E  >  (1) 

m—1  i—1 

where  and  the  expression  for 

i(^)  (u  =  1,  2)  depends  on  the  multi-carrier  scheme  used. 


C.  Channel 

Assume  the  two  antennas  at  both  the  transmitter  and  re¬ 
ceiver  are  sufficiently  spaced,  such  that  the  channels  between 
different  transmit-receive  antenna  pairs  are  independent.  Each 

^Hence,  there  exists  a  one-to-one  correspondence  between  the  pair  (m,  i) 
and  the  sub-carrier  frequency  i . 


TABLE  I 
DEFINITION  OF 


Estimation  phase 


^(fe) 

'^p,q,u,m 


j Su,i,  g  odd 
g  even 


Data  phase 


b(k) 

'^p,q,u,m 


f  f 

Ojp^q^mi  W  ~  O'  odd 

,  tt  =  2,  godd 

—o,p,q,m,  u=  1,  q  even 

«  =  2,  geven 


channel  is  frequency-selective  with  respect  to  the  overall  sys¬ 
tem  bandwidth,  but  each  sub-band  is  assumed  to  be  frequency 
non-selective  with  Rayleigh-distributed  fade  amplitudes.  Con¬ 
sequently,  our  use  of  Alamouti  space-time  block  coding  yields 
fourth-order  spatial  diversity  at  each  sub-carrier  frequency;  the 
theoretical  frequency  diversity  gain,  however,  depends  on  the 
correlation  of  the  channel  gains  across  the  different  sub-bands, 
which  is  a  function  of  the  coherence  bandwidth  of  the  channel. 

In  the  absence  of  Doppler  effects,  there  is  no  inter-channel 
interference  between  adjacent  sub-bands  since  they  do  not 
overlap.  Assuming  perfect  power  control,  the  complex  lowpass 
equivalent  received  signal  at  frequency  fm,i  is  then  given  by 
the  2  X  1  vector  where 


K 


(k) 


s(fe) 


v.u.m.t  u.m. 


k—1  u—1 


(fc) 


v  =  l,2. 

(2) 


v.u.m.i 


A  (fc) 
=  a,'  ' 


,6' 


(fc) 


In  the  above  equation,  g, 
zero-mean  circularly  symmetric  complex  Gaussian  random 
variable  representing  the  fc-th  user’s  channel  gain  at  frequency 
fm  i  between  the  w-th  transmit  antenna  and  the  v-th  receive 

’  (fc') 

antenna,  where  ct),  ^  ^  ^  is  a  Rayleigh  random  variable  with 


unit  second  moment,  and 

,(fc) 


g(fe) 


is  a  uniform  random  variable 


over  [0, 27r).  The  identically  distributed,  but 

their  joint  statistics  depend  on  the  coherence  bandwidth  of 
the  channel.  Also,  represents  the  relative  time  delay 
between  the  /c-th  user  and  the  desired  user,  and  hy{t)  is  a 
zero-mean  complex  additive  white  Gaussian  noise  (AWGN) 
random  process  with  two-sided  power  spectral  density  Nq. 
The  noise  processes  at  the  two  receive  antennas  are  assumed 
to  be  independent. 


D.  Receiver 

The  block  diagram  for  the  receiver  is  shown  in  Fig.  2.  At 
each  antenna,  the  received  signal  first  passes  through  a  bank 
of  MR  demodulators.  We  label  the  output  of  the  demodulator 
at  frequency  i  of  the  u-th  receive  antenna  during  the 
estimation  and  data  phases  as  IT,  „  (<Z  =  ■  j  Qp) 

and  Yq  y  rni  (?  =  Qp  +  l,---)  Q),  respectively.  During  the 

’  ’  ’ 

estimation  phase,  the  estimate  of  ^  j,  Wv^u,m,i,  is  formed 
by  taking  a  sample  average  of  the  {Wq^y^m,i}  during  the  odd 
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Estimation  Phase  (9  =  1,  •  •  • ,  Qp) 


^1,1, 1,1 

1^1, 2, 1,1 
^1,1, M,P 

^2, 1,1,1 
^2, 2, 1,1 


i^2,l,M,P 

^2,2,M,P 


Data  Phase  (g  =  Qp  +  1, . . . ,  Q)  Channel 

Estimates 


Estimated 

Data 


Channel 

Estimates 


Fig.  2.  Receiver  block  diagram. 


or  even  time  slots: 


W, 


v.u.m.t 


=  < 


f 

1 

ilo 

[  1  ] 

9=1 

V  q  odd  j 

1 

ilo 

1  ^  ^ 
gp/2^^™ 

V  q  even  / 

U  =  1, 


u  =  2, 


(3) 


where  the  normalizing  constant  Aq  ensures  that  Wv^u,m,i  is 

(k)  ^  ^ 

an  unbiased  estimate  of  j.  During  the  data  phase,  the 

{Yq,v,m,i}  at  each  sub-carrier  frequency  over  two  consecutive 
time  slots  feed  into  a  linear  Alamouti  decoder,  which  outputs 
the  2x1  vector 

q,m^i  —  -^9+1,771,2]  ,  (4) 

where  q  is  an  odd  integer  in  the  range  Qp  +  1, . . .  ,Q.  In  the 
above  equations,  we  have 


Cr  ™  — 


Wi^l^rn,i 
112,1,777,7 
Wj 


1.2, m,i 


2,2,m,i 


Wi^2, m,i 
112,2,777,7 
-W"i*,l,7„,7 
-11^2*,1,7„,7 


,7  W2,; 


y  q,m,i  — 


^q,l,m,i 

^q,2,m,i 


y 


q-\-l,l,m,i 


q-\-l,2,m,i. 


(5) 


(6) 


such  that  Zq 


and  Z, 


q-\-l,m,i 


We  focus  on  the  detection  of  Zq^rn,i-  After 


Alamouti  decoding,  the  demapper  groups  together  the  R 
variables  associated  with  each  of  the  M  parallel  data  symbols 
and  feeds  them  into  a  bank  of  M  combiners.  The  output  of 
the  TO-th  combiner  is  expressed  as 

I  ^  1 

^q,m  —  7  ^  ^  Yq  jji  i  / 

U=i  J 

=  y  I  l^l,m,^yq,7n,^  +  ^yq,„^,^^l,m,^  \  ,  (7) 

7  =  1  1  1 

where,  again,  q  is  an  odd  integer  in  the  range  Qp  +  1, . . .  ,Q, 
and  the  m-th  data  symbol  during  the  g-th  time  slot  is  decoded 
as  sgn|Zgm|.  Finally,  parallel-to-serial  (P/S)  conversion 
gives  the  estimated  data  sequence. 


E.  Probability  of  Error 
By  defining 

Wm  =  [w^^l  •••  , 

yq,m  =  [5^9,777,1  ■  ■  ■  5^9,777, it] 


(8) 

(9) 


where 


=  [1^1,1,777,7  112,1,777,7  1Pi,2,777,7  112,2,777,7]  )  (10) 

yq^m^i  —  [1)7,1,777,7  1)7,2,777,7  1)7+1,1,777,7  1)7+1,2,777,7]  5  (11) 

we  can  re-express  Zq^ra  in  (7)  as  a  Hermitian  quadratic  form 
in  zero-mean  complex  Gaussian  random  variables. 


Yq^m  -  ''^9,777 777) 


(12) 
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Wi 


I - 1 


W2 


I  I  Symbol  #1  Q  Symbol  #2 


Fig.  3.  Spectra  of  the  MC-DS-CDMA  and  MC-CDMA  signals.  In  this  example,  M  =  2  parallel  data  symbols  are  each  repeated  across  Ri  =  8  and 
R2  =  2  sub-carriers  in  MC-CDMA  and  MC-DS-CDMA,  respectively,  where  TV  =  4  is  the  processing  gain  in  each  MC-DS-CDMA  sub-band.  We  assume 
the  coherence  bandwidth  of  the  channel,  (A/)c,  is  equal  to  IV2. 


where 


Wm 

yq,m 


F 


0477 

J-4R 


O4T? 


(13) 


Assuming  a  -fl  was  transmitted,  the  bit  error  probability  is 
given  by 


P.  =  P- 


^p,q,m  ^  ^ 


,(fc) 


(14) 


A  general  method  for  obtaining  a  closed-form  expression  of 
(14)  is  outlined  in  Appendix  A. 


III.  MC-CDMA  vs.  MC-DS-CDMA 
A.  Overview 

We  now  apply  MC-CDMA  and  MC-DS-CDMA  signaling 
to  the  system  model  presented  in  Section  II.  To  distinguish 
the  system  parameters  between  the  two  schemes,  we  write 
S^Si,W  =  Wi,  and  R  ^  Ri  for  MC-CDMA  and  S  = 
S2,  W  =  W2,  and  R^  R2  for  MC-DS-CDMA.  Each  MC- 
CDMA  sub-band  is  non-spread  and  has  a  symbol  rate  of  l/Tg, 
where  Tg  =  ^MT.  Therefore,  the  bandwidth  of  each  sub¬ 
band  is  Wi  =  l/Tg.  On  the  other  hand,  each  MC-DS-CDMA 
sub-band  has  a  processing  gain  (spreading  factor)  of  N,  so 
W2  =  NWi  =  1/Tc,  where  1/Tc  =  N/Tg  is  the  chip  rate. 
Since  the  two  multi-carrier  schemes  are  constrained  to  the 
same  information  rate  and  bandwidth,  we  have  MRiWi  — 
MR2W2,  which  implies  Ri  =  Wi?2-  An  example  of  how  the 
two  schemes  compare  in  the  frequency  domain  is  given  in 
Fig.  3. 

In  the  following  comparison,  we  assume:  (1)  the  channel  is 
the  same  for  both  schemes,  and  the  coherence  bandwidth  of 
the  channel,  (A/)^,  is  equal  to  W2',  (2)  perfect  carrier,  chip, 
and  symbol  synchronization  are  established  at  the  receiver; 
and  (3)  all  K  users  use  long  spreading  sequences,  which  are 
modeled  as  independent  random  binary  sequences  of  ±l’s. 


B.  MC-CDMA 

1)  Transmitter:  The  block  diagram  of  the  MC-CDMA 
modulator  at  frequency  fm,i  of  the  u-th  transmit  antenna,  as 
shown  in  Fig.  4,  is  identical  to  the  one  in  [20].  Each  symbol  at 
the  input  of  the  MC-CDMA  modulator  is  multiplied  by  only  a 
single  chip  from  the  fc-th  user’s  spreading  sequence,  ^  ^ 
The  resulting  sequence  modulates  an  impulse  train  and  passes 
through  a  chip  wave-shaping  filter  denoted  by  Hi{f  ).  Finally, 


the  output  signal  of  the  filter  modulates  the  corresponding 
sub-carrier. 

The  transmitted  signal  for  the  k-th  user  is  given  by  (1), 
where 


m  =  V  V 

P—  —  00  q—1 

■hi  [t  -  {pQ  -Cq-  l)rs] . 


(15) 


In  the  above  equation,  hi{t)  is  the  impulse  response  of  the  chip 
wave-shaping  filter  iTi(/),  and  Ai^q  is  the  transmit  amplitude 
during  the  g-th  time  slot: 


A-i 


,q 


Ale;  q  —  Qp; 

Aid,  q  =  Qp  -f  1, . . . ,  Q. 


(16) 


To  keep  the  total  energy  transmitted  from  the  two  antennas 
the  same  over  all  time  slots,  we  choose  Aid  =  Aie/-\/2-  As 
a  result,  the  energy  per  bit  can  be  shown  to  equal  Ei,i  = 
\RiA\^  =  RiA^d  [22,  Appendix  P]. 

2)  Channel:  Given  (A/)e  =  W2  =  NWi  »  Wi,  the 
coherence  bandwidth  of  the  channel  spans  multiple  MC- 
CDMA  sub-bands,  which  means  that  the  sub-bands  may  be 
highly  correlated.  We  assume  a  correlated  block  fading  model 
in  which  the  MRi  sub-bands  are  grouped  into  MRi/N 
frequency  blocks — each  block  of  N  sub-bands  corresponds  to 
one  MC-DS-CDMA  sub-band,  and  the  number  of  repetitions- 
per-data  symbol  within  a  block  is  equal  to  i?t,  =  N/M 
(see  Fig.  3).  We  assume  flat  fading  across  each  block,  such 
that  the  N  fade  amplitudes  associated  with  each  block  are 
identical  (i.e.,  perfectly  correlated).  Also,  we  assume  the  fade 
amplitudes  associated  with  any  two  sub-bands  from  different 
blocks  are  independent.  Thus,  even  though  each  symbol  is 
transmitted  across  Ri  sub-carriers,  the  effective  order  of 
frequency  diversity  per  data  symbol  is  only  MRi/N,  and  the 
overall  diversity  gain  per  data  symbol  is  AMRi/N,  where 
the  factor  of  4  comes  from  the  spatial  diversity.  Note  that 
the  frequency  diversity  gain  for  MC-CDMA  does  not  change 
when  we  vary  M,  since  MRi  and  N  are  both  fixed. 

3)  Receiver:  The  complex  equivalent  lowpass  version  of 
the  MC-CDMA  demodulator  in  [20]  is  shown  in  Fig.  4,  where 
hm,i,i{t)  =  h\{—f)  is  the  complex  equivalent  lowpass  impulse 
response  of  the  bandpass  filter.  The  despreading  operation  in 
MC-CDMA  simply  involves  multiplication  by  a  single  chip. 
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Modulator 


Impulse  period  =  Ts 


(j(^) 

p,q,m,i 


Demodulator 

t  =  {pQ  +  q-  l)Ts 


Fig.  4.  MC-CDMA  modulator  at  frequency  fm,i  of  the  ri-th  transmit  antenna,  and  complex  equivalent  lowpass  version  of  the  MC-CDMA  demodulator  at 
frequency  fm,i  of  the  n-th  receive  antenna. 


Modulator 


Impulse  period  =  Tc 


Q 


(fe) 


Demodulator 

t=[{pQ  +  q-l)N  +  n-l]Tc 


Fig.  5.  MC-DS-CDMA  modulator  at  frequency  fm,i  of  tho  u-th  transmit  antenna,  and  complex  equivalent  lowpass  version  of  the  MC-DS-CDMA  demodulator 
at  frequency  fm,i  of  the  f-th  receive  antenna. 


Assume  Hi{f)  is  of  unit  energy,  i.e.,  \Hi{f)\‘^  df  =  1. 

We  define  xi{t)  Xi{f)  =  \Hi{f)\‘^  and  assume 


M.f) 


T _ 3-  <  f  < 

0,  otherwise, 


sin  nt/Tg 
irt/Ts 

(17) 


Since  Xi{f)  satisfies  the  Nyquist  criterion,  inter-symbol  in¬ 
terference  is  not  present. 


C.  MC-DS-CDMA 

1)  Transmitter:  The  block  diagram  of  the  MC-DS-CDMA 
modulator  at  frequency  of  the  u-th  transmit  antenna  is 
identical  to  the  one  in  [20]  (see  Fig.  5).  During  each  time  slot, 
bp,q  u,m  is  multiplied  by  N  chips  in  the  spreading  sequence, 
{ci'^q,n}n=i-  The  transmitted  signal  for  the  fc-th  user  is  also 


given  by  (1),  but  with 

oo  Q  N 

s(^)  (f\  _  \  '  \  '  \  '  A  u{k)  r<{k) 

P—  —  CX)  q—1  n—1 

-  i(pQ  +  q-m  +  n-  1]T,},  (18) 

where  h2  (t)  is  the  impulse  response  of  the  chip  wave-shaping 
filter  H2{f),  and  A2.q  is  the  transmit  amplitude  during  the 
qr-th  time  slot,  which  is  defined  as 


A2,q 


^2e,  9  =  1,  •  ■  •  ,  Qp, 

A2d,  9  =  Qp  +  1,  •  ■  •  ,  Q- 


(19) 


We  choose  A2d  =  A2el'j2,  such  that  the  energy  per  bit  is 
equal  to  E^,2  =  \R2A2eN  =  R2A\^N  [22,  Appendix  P]. 
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2)  Channel:  With  (A/)c  =  W2,  we  assume  flat  fading 
in  each  MC-DS-CDMA  sub-band  and  independent  fading  be¬ 
tween  different  sub-bands.  As  a  result,  the  order  of  frequency 
diversity  seen  by  each  data  symbol  is  i?2,  and  the  overall 
diversity  gain  per  data  symbol  is  4i?2-  Since  i?i  =  NR2,  we 
have  4:MRi/N  =  4Mi?2  >  4i?2,  which  implies  that  MC- 
CDMA  has  M  times  the  diversity  of  MC-DS-CDMA. 

3)  Receiver:  We  use  the  complex  equivalent  lowpass  ver¬ 
sion  of  the  MC-DS-CDMA  demodulator  in  [20],  where 

(see  Fig.  5).  In  contrast  to  MC-CDMA, 
the  despreading  operation  in  MC-DS-CDMA  involves  multi¬ 
plication  by  N  chips.  We  define  X2{t)  X2{f)  =  |iF2(/)P, 

where  the  expression  for  X2(f)  is  given  by  replacing  Tg  with 
Tc  in  (17).  Since  X2{f)  also  satisfies  the  Nyquist  criterion, 
inter-symbol  interference  is  not  present. 

IV.  Performance  Analysis 
In  this  section,  we  present  the  expressions  for 
Yq,v,m,i,  and  Wv^u,m,i  as  they  apply  to  MC-CDMA  and  MC- 
DS-CDMA.  We  also  give  closed-form  results  for  the  bit  error 
probabilities  of  both  schemes,  assuming  the  desired  user  is 
k  =  1  (hence,  =  0). 


A.  MC-CDMA 


1 )  Demodulator  Output:  The  output  of  the  demodulator  at 
frequency  fm,i  of  the  u-th  antenna  during  the  g-th  time  slot 
is  given  by 


—  Ale 

g  =  1, . . . ,  Qp, 

2 

—  Aid  E  ^p,q,u,m  9v,u,m,i 

u—1 

q  =  Qp  +  1,  •  ■  • ,  Q, 


(20) 

(21) 


where 


Iw„. 


=  a 


K  2 

k^2  u^l 
00  Q 


/  j  /  j  p' ,q' ,u^m^p' ,q' ,m4 

p'  —  —  oo  q' —  1 

■  Xi[t-  (p'Q  +  q'  -  l)Ts 


r(fe)l 


t={pQ+q-l)Tg 


(22) 


and 


t=(pQ+q-l)Tg 


(23) 

represent  the  multiple  access  interference  (MAI)  due  to  the 
other  K  —  1  users  and  the  contribution  of  the  AWGN  to  the 
test  statistic,  respectively.  The  expressions  for 

n.  i 


are  identical  to  (22)  and  (23),  respectively,  except 


they  are  associated  with  a  different  range  of  values  for  q.  Note 
that  the  desired  signal  and  AWGN  terms  are  both  zero-mean 
complex  Gaussian.  Thus,  in  the  single-user  case,  Wg  „  and 
are  both  zero-mean  complex  Gaussian  random  vari¬ 
ables.  When  K  is  large,  the  MAI  terms  are  approximated  as 


zero-mean  complex  Gaussian  random  variables,  and 

and  are  approximated  as  zero-mean  complex  Gaussian 

random  variables  as  well. 

2)  Channel  Estimation:  Substituting  (20)  into  (3)  and  set¬ 
ting  Aq  =  Ale,  the  estimate  of 


W, 


v,u,m,z 


9 


(1) 

v,u,m,i 


< 


9 


(1) 

v,u,m,i 


1  1 

Ale  Qp/2 

1  1 

Ale  Qp/2 


Qp 

E 


9=1 

q  odd 

Qp 


E 


,m,i  ’ 


9=1 

q  even 


tt  =  1, 

U  =  2, 
(24) 


where  Twg,g,„,,i  =  Iwg.v.m.i  +  Since  Wv^u,m,i  is 

a  linear  combination  of  the  it  is  also  a  zero- 

mean  complex  Gaussian  random  variable  when  K  =  1  and 
is  approximated  as  a  zero-mean  complex  Gaussian  random 
variable  when  K  is  large. 

From  estimation  theory,  both  {Wq^y^m,i}q  oM  and 
{Wq,v,m,i}qeYsn  represent  linear  models.  The  classical 
approach  assumes  j  is  deterministic,  and  we  can  show 

that  is  the  minimum  variance  unbiased  (MVU) 

estimator  [23].  On  the  other  hand,  the  Bayesian  approach 
assumes  ^  ^  is  zero-mean  complex  Gaussian.  It  can  be 
shown  that  the  minimum  mean  square  error  (MMSE)  estimate 
of  ^  ^  is  just  a  scaled  version  of  Wy^u,m,i,  and  that  both 
the  MVU  and  MMSE  estimators  yield  the  same  probability 
of  error  [22,  Appendix  E]. 

3)  Probability  of  Error:  Using  the  approach  outlined  in 
Appendix  A,  the  probability  of  error  is  derived  as  (see 
Appendix  B) 


Pe  = 


b  b-l 

E 

i2—0 

1  —  931 


6-1-1-  72^  f  1 
*2 


-f 


1  -  9A1 
2 

1  —  913 
2 

1  —  923 
2 

1  —  943 


a  6—1  —  72 

E 

*3=0 
a  b— 1  —  12— is 


E 

r- 
-1 

E 

7l=0 
h a— 1—7 

E 


r—Q 
b  a— 1 


0-1-1-  73^  1  +  931 

*3 

a— l  +  r^  f  I  +  941 
r 


6-1-1-  7i\  1  +  913 

ii 


72=0 


1  ^6-1+72\  (1+923'^"^ 
*2 


a  a— 1  —  71—72 


E 


r=0 


a  —  1  -\-r 
r 


1  +  943 


(25) 


where  a  =  4Mi?i/At,  6  =  4(i?i  -  MRi/N), 


Rh  — 


931  = 


Rb  ^Jd(^  Yo  )  +  (  No  ) 


(26) 
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M41  = 


M13  = 


M23  = 


M43  = 


Rb  +  ^ 

ld\ 

Rb  +  ^ 

fd\ 

(S) 

1  - 1 

ld\ 

(tt) 

0 

1  + 1 

[d\ 

(^) 

0 

[%] 

1  +  1 

ld\ 

(^) 

0 

[d\ 

(^) 

0 

Rb 


<t-) 


(27) 


(28) 


(29) 


(30) 


In  the  above  equations, 


we 


+  1  “1“ 


V 


Qp  1^ 


Cl 


(^) 


have 
■f  2 


Qp  1^ 


Cl 


(^)r. 


where  ci(^j  =  (K  -  1)  + 

4)  Asymptotic  Analysis — Perfect  CSI:  While  the 

complexity  of  the  above  expressions  makes  it  difficult 
to  gain  an  intuitive  understanding  of  the  probability  of 
error,  the  usefulness  of  this  closed-form  result  becomes 
apparent  when  we  consider  the  asymptotic  case  of  perfect 
CSI.  Suppose  we  have  a  completely  static  channel,  such  that 
the  fade  amplitudes  are  constant  over  an  infinitely  long  data 
frame  (i.e.,  Q  oo).  If  we  let  Qp  oo,  then  the  error 
terms  in  (24)  vanish,  and  the  channel  estimation  becomes 
perfect.  Since  limQp_>oo 
we  have  limQp_>oo  Msi 
liruQp—^oo  /ri3  ~ 
liruQp—^oo  /r43  — 


(  tVo  )  “  (  JVoO  )  ’ 

=  limQp_>oo /r4i  =  1, 


lirnQp_»oo  p,23 
Rh _ 


7p  = 


_  Rb  _ 


Pb+2ci( 
1 


^bi 
^0  , 


=  —1,  and 


where 


1+7P  ’ 


No 


1  f  1  Pm  ^ 

2  \MRi/N  No  ) 


these  results,  we  finally  obtain 


lim  Pe  = 

Qp^oo 


1-^ 


a  a—1 


2=0 


a—1 

i 


1  +  /X 


Using 


(31) 


where  P  =  y  1+7'  (^0  is  simply  the  probability 

of  error  of  an  ideal  maximal  ratio  combiner  with  AMRi/N- 
th  order  diversity  and  an  average  signal-to-interference-plus- 
noise  ratio  (SINR)  per  diversity  branch  of  fp  [24,  Eqs.  (14.4- 
15)-(14.4-16)]. 


B.  MC-DS-CDMA 

1)  Demodulator  Output:  The  demodulator  output  at  fre¬ 
quency  fm,i  of  the  u-th  antenna  during  the  g-th  time  slot  is 
given  by 

q  =  l,...,Qp,  (32) 

2 

Yq,v,m,i  =  ^2dA/Vy^  ^p}q,u,m  9v}i,m,,i  3"  N.u.m.i 
u—1 

<Z  =  Qp  +  l,...,g,  (33) 


where 


K  2 


1  iv  IS.  ^  oo  Q 

=  -7=  ^  ^  9^vl,m,^  XI  X 

^  n—1  k—2u—l  p'  =  —  ooq'  —  l 


and 


N 


fc=2  ii=l 
L  (fc)  /^(^) 


/  j  ^p' ,q'  ,u,m^p'  ,q' ,n' 

n'  —  l 

■X2{t  -  [{p'Q  -\-q'  -  1)N  -\-n'  -  l]Tc 


-r«} 


*=[(pQ-|-<?— 1)^  +  "— llT: 


(34) 


N 


n=l 

•  (ny{t)  *  hm,i,2 


t=[{pQ-i-q—l)N +n—l]Tt, 


(35) 


correspond  to  the  MAI  and  AWGN  terms,  respectively.  Again, 
the  expressions  for  are  identical  to 

(34)  and  (35),  respectively.  As  in  the  MC-CDMA  analysis, 
Wq^y^jn,i  and  Uj  t,  m  i  are  both  zero-mean  complex  Gaussian 
random  variables  in  the  single-user  case  and  are  approximated 
as  zero-mean  complex  Gaussian  random  variables  when  K  is 
large. 

2)  Channel  Estimation:  The  expression  for  is 

obtained  by  substituting  (32)  into  (3).  For  MC-DS-CDMA, 
we  require  a  normalizing  constant  of  Aq  =  A2e'/N- 

3)  Probability  of  Error:  From  Appendix  C,  the  probability 
of  error  is  derived  as 


where 


(37) 


(38) 


and  Note  that  this  is  just  the 

probability  of  error  of  an  ideal  maximal  ratio  combiner  with 
4i?2-th  order  diversity,  except  that  the  average  SINR  per 
diversity  branch  has  been  replaced  by  7  of  (38)  [24,  Fqs.  (14.4- 
15),(14.4-16)]. 

4)  Asymptotic  Analysis — Perfect  CSI:  Once  again,  the  re¬ 
sults  for  perfect  CSI  can  be  obtained  by  assuming  both  a 
static  channel  and  an  infinitely  long  data  frame  {Q  00) 
and  then  computing  the  limits  of  the  expressions  in  (36)- 
(38)  as  Qp  — ^  00.  Since  Qp  is  only  present  in  7,  we  can 
simply  replace  7  in  /x  with  fp  =  limgp^oo  7  =  — = 

1  /  1 

2  \R2  No  j  1  I  K-i  1 

\  y  7J-  R2  JVq 
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V.  Numerical  Results 

A.  Review  of  General  Trends 

We  briefly  summarize  the  general  trends  observed  in  Sce¬ 
nario  #2  of  [20]  for  the  SISO  system.  When  M  >  1  and 
K  —  1,  the  error  probability  curves  for  MC-CDMA  and  MC- 
DS-CDMA  always  cross.  This  cross-over  point,  which  we  de¬ 
note  by  (A,Pe,A),  implies  that  Pe.mc-ds-cdma  <  Tb.mc-cdma 
for  Eb/No  <  A  and  Pe.mc-ds-cdma  >  Tb.mc-cdma  for 
Eb/No  >  A.  Likewise,  we  also  have  (P6/Wo)mc-ds-cdma  < 
(Pb/Ao) MC-CDMA  for  Pe  >  Pe,A  and  (Pb/Ao)MC-DS-CDMA  > 
(Pb/Ao)MC-CDMA  for  Pe  <  Pe,A-  When  M  >  1  and  a:  >  1, 
the  MC-DS-CDMA  and  MC-CDMA  curves  cross  only  if  Qp 
is  sufficiently  large;  otherwise,  MC-DS-CDMA  gives  a  lower 
Pe  for  all  Eb/Nf).  This  is  due  to  the  fact  that  our  detection 
schemes  only  involve  matched  filtering  and  lack  a  mechanism 
for  combatting  MAI.  As  a  result,  an  error  floor  is  present  at 
high  Eb/Nf),  and  the  curves  for  both  multi-carrier  schemes 
flatten  out. 

In  the  figures  to  be  presented  shortly,  we  will  see  that  these 
trends  for  the  SISO  system  also  apply  to  our  MIMO  system. 
This  is  because  the  multi-carrier  comparisons  in  both  cases 
are  centered  around  the  trade-off  between  diversity  gain  and 
channel  estimation  errors. 

B.  SISO  Vi.  MIMO  Comparison 

Using  the  closed-form  expressions  for  the  bit  error  proba¬ 
bilities  of  MC-CDMA  and  MC-DS-CDMA  presented  in  the 
previous  section  and  in  Section  III-E  of  [20],  we  compare  the 
MIMO  results  with  those  of  Scenario  #2  of  the  SISO  system 
to  study  the  impact  of  the  additional  diversity  on  the  multi¬ 
carrier  comparison.  For  both  the  SISO  and  MIMO  systems,  we 
assume  Si  =  MRi  =  512  and  S2  =  MR2  =  8  sub-carriers 
for  MC-CDMA  and  MC-DS-CDMA,  respectively,  where  the 
processing  gain  of  each  MC-DS-CDMA  sub-band  is  A  =  64. 
We  set  Ebi  =  Eb2  =  Eb  and  consider  both  K  =  1  and 
K  =  16.  The  overall  diversity  gains  per  data  symbol  in  MC- 
DS-CDMA  and  MC-CDMA  are  equal  to  R2  and  MRi/N, 
respectively,  in  the  SISO  system,  and  4i?2  and  AMRi/N, 
respectively,  in  the  MIMO  system.  As  we  pointed  out  earlier 
in  Section  III,  MC-CDMA  has  M  times  the  diversity  of  MC- 
DS-CDMA.  Thus,  when  M  =  1,  both  multi-carrier  schemes 
have  equal  diversity,  and  the  probability  of  error  was  shown 
to  be  always  lower  in  MC-DS-CDMA  due  to  more  reliable 
channel  estimation  [20].  Henceforth,  we  only  consider  M  >  1. 
Since  the  transmission  of  pilot  symbols  during  the  estimation 
phase  occurs  during  every  time  slot  in  the  SISO  system,  but 
only  during  every  other  time  slot  at  each  antenna  in  the 
MIMO  system,  we  assume  Qp  =  16  for  the  SISO  system 
and  consider  both  Qp  =  16  and  Qp  =  32  for  the  MIMO 
system.  When  Qp  =  16,  the  MIMO  system  maintains  the 
same  throughput  as  the  SISO  system®  (i.e.,  both  systems 
transmit  Qp  data  symbols  per  frame  at  each  sub-carrier),  but 
when  Qp  =  32,  the  MIMO  system  uses  the  same  amount  of 
total  energy  (i.e.,  same  number  of  pilot  symbols)  as  the  SISO 
system  to  estimate  each  channel  gain. 

assume  the  SISO  and  MIMO  systems  use  data  frames  of  the  same 
length  (i.e.,  same  Q). 


Fig.  6.  SISO  vs.  MIMO:  K  =  1,  M  =  2,  R2  =  4:,  MRi/N  =  8.  The 
orders  of  diversity  (per  data  symbol)  in  MC-DS-CDMA  and  MC-CDMA  are 
equal  to  4  and  8,  respectively,  in  the  SISO  system,  and  16  and  32,  respectively, 
in  the  MIMO  system. 

First,  we  consider  the  single-user  scenario  and  compare  the 
SISO  and  MIMO  results  against  a  target  bit  error  rate  of 
P*  =  10“®.  The  case  M  =  2  is  examined  in  Fig.  6.  Here,  MC- 
CDMA  has  twice  the  frequency  diversity  of  MC-DS-CDMA. 
Note  that  the  MC-DS-CDMA  and  MC-CDMA  curves  for  the 
SISO  system  cross  near  P*,  which  means  that  both  multi¬ 
carrier  schemes  require  approximately  the  same  Eb/Ng  to 
achieve  P*.  For  the  MIMO  system,  however,  we  observe  that 
for  both  Qp  =  16  and  Qp  =  32,  the  Eb/No  values  required 
to  achieve  P*  are  lower  in  MC-DS-CDMA  by  at  least  3dB. 
Even  if  we  reduce  the  target  bit  error  rate  to  P*  =  10“^, 
the  required  Eb/No  values  are  still  lower  in  MC-DS-CDMA. 
Next,  we  consider  M  =  8.  We  observe  in  Fig.  7  that  the 
cross-over  points  in  the  SISO  case  and  in  both  instances  of 
the  MIMO  case  occur  at  Pe  values  which  are  greater  than 
P*  =  10-5.  Thus,  MC-CDMA  achieves  P*  at  lower  Eb/No 
values  than  MC-DS-CDMA  in  both  the  SISO  and  MIMO 
cases,  although  the  differences  in  the  required  Eb/No  values 
are  much  smaller  in  the  two  MIMO  cases.  Based  on  these 
results  for  M  =  2  and  M  =  8,  we  conclude  that  for  K  =  1, 
the  performance  gains  obtained  by  MC-DS-CDMA  from  the 
addition  of  spatial  diversity  are  enough  to  offset  the  difference 
in  frequency  diversity  between  the  two  schemes  when  M  is 
small,  but  not  when  M  is  large. 

Now,  consider  a  multi-user  scenario  with  K  =  16.  When 
M  =  2,  we  see  in  Fig.  8  that  the  MC-CDMA  and  MC-DS- 
CDMA  curves  do  not  cross  in  either  the  SISO  or  MIMO  cases. 
Furthermore,  at  a  given  Eb/No,  the  differences  in  probability 
of  error  between  the  two  schemes  are  much  larger  in  the  two 
MIMO  cases  than  in  the  SISO  case.  This  implies  that  in  order 
for  the  MC-CDMA  and  MC-DS-CDMA  curves  to  cross,  we 
will  have  to  use  a  larger  Qp  in  the  MIMO  case  than  in  the 
SISO  case.  When  M  =  8,  the  curves  for  the  two  schemes  in 
Fig.  9  do  cross  in  the  case  of  the  SISO  system  when  Qp  =  16. 
In  the  MIMO  system,  however,  a  cross-over  point  does  not 
exist  for  either  Qp  =  16  or  Qp  =  32;  we  must  use  Qp  >  32 
in  order  for  the  curves  to  cross.  Thus,  we  conclude  that  the 
additional  diversity  increases  the  value  of  Qp  required  for  the 
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Fig.  7.  SISO  vs.  MIMO:  K  =  1,  M  =  8,  R2  =  I,  MRi/N  =  8.  The 
orders  of  diversity  (per  data  symbol)  in  MC-DS-CDMA  and  MC-CDMA  are 
equal  to  1  and  8,  respectively,  in  the  SISO  system,  and  4  and  32,  respectively, 
in  the  MIMO  system. 


Fig.  8.  SISO  vs.  MIMO:  K  =  16,  M  =  2,  R2  =  4,  MRi/N  =  8. 


MC-CDMA  and  MC-DS-CDMA  curves  to  cross  and  for  there 
to  be  a  trade-off  in  performance  between  the  two  multi-carrier 
schemes. 


Fig.  9.  SISO  vs.  MIMO:  K  =  16,  M  =  8,  R2  =  1,  MRi/N  =  8. 


offset  the  difference  in  frequency  diversity  between  the  two 
schemes  when  M  is  small,  but  not  when  M  is  large,  where 
M  represents  the  number  of  parallel  data  streams.  Also,  when 
matched-filter-based  detection  is  used  for  the  case  of  multiple 
users,  our  results  showed  that  the  additional  diversity  increases 
the  number  of  pilot  symbols  required  to  force  a  performance 
trade-off  between  the  two  schemes  (i.e.,  to  get  the  MC-CDMA 
and  MC-DS-CDMA  curves  to  cross). 

Appendix  A 

Hermitian  Quadratic  Forms  in  Complex  Gaussian 
Random  Variables 

Consider  the  general  Hermitian  quadratic  form 

Z  =  v^Fv,  (39) 

where  v  =  [ui  •  •  •  is  a  column  vector  of  N  jointly 

distributed  zero-mean  complex  Gaussian  random  variables, 
and  F  is  an  A^  X  Af  Hermitian  matrix.  Define  the  N  x  N 
covariance  matrix  of  v  as 

R  =  f  f;  [vv^]  .  (40) 


VI.  Conclusion 

Our  previous  comparison  between  MC-DS-CDMA  and 
MC-CDMA  focused  on  the  trade-off  between  diversity  gain 
and  channel  estimation  errors,  but  only  a  SISO  system  was 
considered.  In  this  work,  we  incorporated  spatial  diversity  into 
the  comparison  by  using  a  MIMO  system  employing  Alamouti 
space-time  block  coding.  To  quantify  the  effects  of  this 
additional  diversity,  we  derived  closed-form  expressions  for 
the  bit  error  probabilities  of  the  two  multi-carrier  schemes,  and 
we  compared  the  results  of  the  MIMO  system  against  those 
of  the  SISO  system.  Since  we  only  considered  those  cases 
where  MC-CDMA  has  higher  frequency  diversity  than  MC- 
DS-CDMA,  we  argued  that  the  additional  diversity  benefits 
MC-DS-CDMA  more  than  MC-CDMA.  It  was  shown  for  the 
case  of  a  single  user  that  these  gains  for  MC-DS-CDMA  can 


The  characteristic  function  of  Z,  $2  (iv)  =  E\e^^^],  is  given 
by  [25,  Eq.  (B-3-16)] 


^z{jv) 


1 
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(1  -  2juA„)“™ 

n—1 


A  rrin 

EE 


n=l 


fc+1) 

(1  -  2ji;A„)'=  ’ 


(41) 


where  Ai , . . . ,  Aa  are  the  A  <  N  —  toq  distinct  non-zero 
eigenvalues  of  the  matrix  RF  with  multiplicities  mi, . . . ,  mA 
such  that  iTii  =  Af  —  mo,  and  toq  is  the  number  of 

eigenvalues  of  RF  equal  to  zero.^  The  last  equality  in  (41)  is 
obtained  by  applying  a  partial  fraction  expansion  (PEE).  The 


^The  eigenvalues  equal  to  zero  only  contribute  a  factor  of  1  to  'l’z(fu). 
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PFE  coefficient  Cn^k  can  be  written  in  closed-form  as  [26] 
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As  a  result,  we  have  the  partitioned  matrix 


RF  = 


1 


S 

S 


ww 

H 

wy. 


and  its  characteristic  polynomial  can  be  expanded  as  [23, 
Sec.  Al.1.3] 


(42) 
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where  Qn  denotes  the  set  of  integers 
{ii,  ■  •  ■  5  ^n— 15  ^n+ii  ■  •  ■  5  ^a}  such  that  ip  —  k  1. 

By  defining 
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We  show  in  [22,  Appendix  M]  that  the  matrices  'Sww, 


and  are  derived  as 
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Assuming  binary  signaling,  the  probability  of  error  is  derived 
as  [27,  Eqs.  (46)-(48)] 

^0  1  ^00 
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where  is  the  Ry  x  Rb  all-ones  matrix,  ci  =  {E ' 
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(45) 


In  general,  both  the  eigenvalues  of  RF  and  the  {c„  fc}  will 
need  to  be  evaluated  numerically.  But  if  A  is  small,  then  (44) 
may  reduce  to  a  manageable  closed-form  expression,  and  it 
may  also  be  possible  to  derive  the  eigenvalues  of  RF  in 
closed  form.  For  example,  when  A  =  2,  the  probability  of 
error  reduces  to  [24,  Eq.  (14.4-15)]. 

Appendix  B 

Derivation  of  Probability  of  Error;  MC-CDMA 

We  derive  the  probability  of  error  for  MC-CDMA  using  the 
approach  given  in  Appendix  A.  We  assume  a  4-1  is  transmitted 
(i.e.,  a'p}q,m  =  +1),  and  we  condition  on  From  (40) 

and  (13),  the  covariance  matrix  of  „  is  equal  to 


Substituting  these  expressions  into  (51),  we  can  show  after 
much  algebraic  manipulation  that  the  characteristic  polyno¬ 
mial  of  RF  reduces  to  [22,  Appendix  M] 

det(Al8fli  -  RF)  =  ■  [T{X)f^E/N  ^ 

(56) 

where 


W{X)  =  A^  -  ^ 


8i?i  Qp/2 


Cl 


\No 

T{X)  =  X*  +  bX^  +  cX^  +dX  +  e, 


b  =  -J^Rb, 

Ki 


(57) 

(58) 

(59) 


c  = 


(46) 

(47) 

(48) 

(49) 

(50) 


Ebl  I 

4i?i\ 


R?  —  Rb  (  1  + 


QpI‘2‘ 


Cl 


£^61 


QpI‘2- 


(Ebl 

''\No 


I  Ebl  Ebl  n  f  Ebl 
1 


1  + 


QpI‘2 


1 

Qp/"2 


(60) 


(61) 


_(EbiV\  (Ebi^^'" 
®  (sRi)  ""HWo 


1  2 


■(62) 


Authorized  licensed  use  limited  to:  Univ  of  Calif  San  Diego.  Downloaded  on  February  4,  2010  at  1 1 :00  from  IEEE  Xplore.  Restrictions  apply. 


LING  and  MILSTEIN:  THE  EFFECTS  OF  SPATIAL  DIVERSITY  AND  IMPERFECT  CHANNEL  ESTIMATION  ON  WIDEBAND  MC-DS-CDMA  AND  MC-CDMA  2999 


Since  the  eigenvalues  of  RF  are  given  by  the  solutions 
of  det(Al8fli  —  RF)  =  0,  we  first  solve  for  the  roots  of 
along  with  their  corresponding  multi¬ 
plicities,  to  obtain 
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Likewise,  using  the  general  method  for  solving  quartic  poly¬ 
nomials  [28,  Sec.  2.3],  the  solutions  to  =  0  are 

given  by 


(67) 


(68) 


and 


where  d 


m3  =  AMRi/N, 
TO4  =  AMRi/N, 


i^) 


=  + 


2Rb(^ 


1-f 


1 

Qpf^ , 


o  1 
^'QEn 


Cl 


where 


No  J  N 


R2 


and  B2(e)  is  defined  in 
(55).  Subkitufing  these  expressions  into  (51)  and  replacing 
Ri  with  R2,  we  can  show  that  the  characteristic  polynomial 
of  RF  reduces  to  [22,  Appendix  L] 


det(Al8fl2  -  RF)  =  (A^  -bX-c) 


4it2 


where 


(75) 


(76) 


(77) 


The  equation  det(Al8_R2  —  RF)  =  0  yields  A  =  2  distinct 
eigenvalues,  both  of  multiplicity  mi  =  m2  =  4i?2: 


Ai  — 


b  —  'Jb'^  -f  4c 


A.  =  (78) 


For  >  0,  we  have  C2  which  implies  6^  -|-4c  > 

b^.  Substituting  (44)  into  (45)  and  using  the  fact  that  Ai  is  the 
only  negative  eigenvalue,  we  can  finally  show  that 


(69) 

4R2 

(70) 

^^  =  E 

)  + 

k^l 

/  1 

.  Thus,  RF  has  A  =  4  distinct  eigenval¬ 
ues.  For  >  0,  we  have  ci^-^^  >  0,  which  implies 

>  R^.  Ai  and  A3  are  then  the  only  negative  eigen¬ 
values,  and  the  expression  for  the  probability  of  error  in  (45) 
reduces  to 

4(iti-Mi?i/Ar)  4MRi/N 

Pe=  X!  +  X!  C3,fc.  (71) 


Cl,fc 


1  —  /421 


4^2  4it2-l 

E 

1  +  M2I 


where  /121  =  = 


-\/b^+4c 

shown  to  reduce  to  the  expression  in  (38) 
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4c 


(79) 


can  be 


By  substituting  (44)  into  (71),  we  can  show  after  a  series  of 
manipulations  that  is  equal  to  (25)  [22,  Appendix  M].  The 
expressions  for  p2i,  M3i>  M4i>  Mi3>  d23,  and  ^43  are  obtained 
by  using  (63)-(64)  and  (67)-(68)  in  (43),  and  we  can  easily 
verify  that  ^21  =  0. 


Appendix  C 

Derivation  of  Probability  of  Error:  MC-DS-CDMA 


The  probability  of  error  for  MC-DS-CDMA  is  also  derived 
using  the  approach  presented  in  Appendix  A.  Again,  we 
assume  =  +1,  and  we  condition  on  OpEi  rn-  The 

expressions  for  R,  ^yy,  ^wy,  and  RF  in  (46)-(50) 

still  hold.  For  MC-DS-CDMA,  the  matrices  '^yy,  and 

are  derived  as  [22,  Appendix  L] 
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